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ON LOCAL NON-ZERO CONSTRAINTS IN PDE WITH 
ANALYTIC COEFFICIENTS 

GIOVANNI S. ALBERTI AND YVES CAPDEBOSCQ 


Abstract. We consider the Helmholtz equation with real analytic coefficients 
on a bounded domain Q C We take d + 1 prescribed boundary conditions 
and frequencies a; in a fixed interval [a,b]. We consider a constraint on 
the solutions of the form ^ • • • ? ^ 0, where ^ is 

analytic, which is satisfied in Q. when oj = 0. We show that for any Q' 
and almost any d + 1 frequencies in [a, 6], there exist d + 1 subdomains 
Qk such that O' C UfcOfc and '''' ’ 7^ 0 in Ofe. 

This question comes from hybrid imaging inverse problems. The method used 
is not specific to the Helmholtz model and can be applied to other frequency 
dependent problems. 


1. Introduction 


The motivation for this work comes from hybrid, or multi-physics, parameter 
identification problem in boundary value problems for partial differential equations 
mii]- In such imaging modalities, one part of the inverse problem can be described 
in general terms as follows. Suppose that are the solutions of a 

partial differential equation of the form 

f P{x,u},ul) =0 in O, 

1 < = r on dn, 


for i = 1,..., N. Suppose further that the parameter uj is known, /* is known, 
and some pointwise information is known about a functional of the solutions, e.g. 

..., ..., Wu^), in 11', a subdomain of If, or possibly in all of fl. 

The problem is to reconstruct spatial dependence of the operator P from this in¬ 
formation. In this article we focus on a particular model, a problem of Helmholtz 
type, in a smooth bounded domain 17 C d > 2, given by 


J —div(a Vu) — -I-iwtj) u = 0 in 17, 

\ u = f on (917. 


We assume that a G L°° (I7; and that a is symmetric and uniformly positive 

definite and bounded, that is, for all ^ S there holds 

(2) 1^1^ < < A 1^1^ a.e. in 17 


for some positive constant A, whereas e,a G L°°(I7;M) satisfy 
(3) A“^ <e<A, 0<cr<A a.e. in 17. 
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Assumptions (i) and f guarantee that problem 0 has a unique solution in 
(ri; C) for every f G (9^1; K) and a; G Z? = C \ S, where E denotes the set of 
the discrete Dirichlet eigenvalues of the problem. Let A = [Amin, A^ax] represent 
the frequencies we have access to, for some 0 < Amin < Amax- For simplicity, we 
suppose that A C £). 

Given several boundary conditions /* G and frequencies oJk G A, the 

pointwise information available (or observable) in this case could be 

ul,, or or 

for some i,j and k. As all these data come from measurements, it is important 
to know a priori that the numbers obtained are not mostly measurement error or 
background noise: we want to ensure that the modulus of these data is non-zero. 
For instance, given w and a boundary condition /^, we want to ensure that 

(4) ui ^ 0. 

Alternatively, combining multiple data and constraints into one functional, given w 
and d -I- 1 boundary conditions /^,..., we write 

(5) J^(a;,/\...,/‘^+^) =uidet 

and we want to ensure that 

(6) in A'. 

The constraints (|^ and (|^, or related quantities, appear in laiiniiiiiiijE]. If 
w is large, that is, greater than a constant depending on A and fl only, Q (and a 
fortiori 0) cannot be satisfied in the whole domain, as for any boundary condition 
in the field must cancel (at least when cr = 0). Thus multiple 

boundary conditions or frequencies must be considered. 

Definition 1. Given a finite set of frequencies {uji,... ,ujk} G A^ and a finite 
set of boundary conditions /^,..., (90; K), we say that {wi,..., ujk} x 

{/^,..., /^} is a set of measurements. 

The first concern is whether there exists a set of measurements such that ([^ is 
satisfied everywhere by a subset of this set. The precise meaning of this statement 
is given by the following definition. 

Definition 2. Take fl' C fl. Given K,N G N*, a set of measurements {uJi ,..., Wic}x 
{/^,..., f^} is F-complete in il' if there exists an open cover of fl' 

= u 

P=1 

such that for each p there exist k G {1,..., AT} and ii,..., G {1,..., ZV} such 
that 

(7) \F{u;k,rF..,f^+^){x)\>0, x&%. 

In other words, a A-complete set of measurements gives a cover of the domain 
n' into a finite collection of subdomains, such that the constraints 0 are satisfied 
in each subdomain for different frequencies and boundary conditions. 




VA 


Vu,'^+i 
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Several results show the existence of such i^-complete sets. In the single-frequency 
case, namely with K = \, the existence of such F-complete sets can be proved by 
using Complex Geometric Optics (CGO) solutions |16lE01ITB] or the Runge approx¬ 
imation property HHHHH] under appropriate regularity hypotheses on a, e and 
a. When using CGO, only P = 2 subdomains are needed, while with the Runge 
approximation approach P is larger than 2. These approaches do not indicate how 
suitable boundary conditions should be chosen in practice, as the proof of their 
existence relies on the unknown coefficients: without additional a priori informa¬ 
tion, the search for these boundary conditions requires many trials. Moreover, since 
CGO are very oscillatory, they may require a very sophisticated practical apparatus 
to be implemented. 

An alternative method consists in fixing the boundary conditions and varying 
the frequency. 


Theorem 3. Take d = 2 and suppose that fl is convex and a G 0°’“ for some 
a > 0. There exist K G N* and C > Q depending only on O, O', A, a, Amin ond 
Amax such that 

k — 1 

~ Amin T (Amax Amin) TY T : A: = 1, . . . , AT} X {1, Xi, X 2 } 

A — i 

is F-complete in O'. More precisely, there exists an open cover O' = OJ, such 
that 

\F{uJk, 1, xi,X 2 )\ > C in Tt'f.. 

This result also holds when d = 3 with the boundary conditions l,xi,X 2 ,X 3 provided 
that |ja — Id\\(^o,a < S where 6 depends only on O, A and a. 


This result is proved in j21 (see also mm). In this result, K is bounded a 
priori and possibly large, but the boundary conditions are fixed a priori and non- 
oscillatory. The proof of Theorem relies on an analytic continuation argument 
with respect to the frequency, immersed in the complex plane. When w = 0, 
problem Q becomes 

J —div(a Vuo) =0 in O, 

\ uq = f on 90. 

For the boundary condition /^ = 1, the solution is simply Uq = 1. In two dimen¬ 
sions, for the boundary conditions = Xi and = X 2 , since O is convex, it is 
known that det (Vug, Vug) > 0 in O musiis]. Therefore 


( 8 ) 


F(0,1, xi,X 2 ) = 1 X det 


1 ul ul 

0 VUg 


> 0 in 0. 


In higher dimension such a result is not available, however Schauder elliptic regu¬ 
larity theory shows that if 


(9) ||a — IdWcO’^ < 

for some d > 0 small enough, then 

■ 1 ul ... ug+i 

Vug 


(10) ^(0,1, xi,..., Xd) = 1 X det 


0 


> 0 in O. 
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The principle of the proof of Theorem|^is then to show that this positivity property 
can be transported to any interval ^ in a predictable manner |2] • 

The goal of this article is to investigate what is the minimal number of required 
frequencies K (or, equivalently, the number P of subdomains in Definition in 
dimension d > 2, for the fixed boundary conditions l,Xi,... ,Xci. We consider a 
technically very convenient particular case, namely we assume that 

(11) a,e, and tr are real analytic in fl. 

The main result of the paper reads as follows. 


Theorem 4. Assume that § and ( |11[ ) hold true. Suppose that A Q D and 
take F as in ([^, fl' (s and /^,..., G Hi (5r2;M). If 


(12) F{0j\...,f+^){x)^0, x€n, 

then 


{(wi,... ,a;d+i) e : {uJk}k x {/\ ■ ■ ■,is F-complete in D'} 
is open and dense in . 


Remark 5. This argument can be used for many other partial differential equations, 
see |2] and umiiais] for variants of this argument. It applies in particular to the 
anisotropic Maxwell system of equations. 


Remark 6. There is nothing special about the function F, which was used as an 
illustration in this paper. We only use that F{uj, 1, cci,..., Xd) expressed in terms 
of ul,.. .,ui+\ that is, 

F{ujA,Xi,...,Xd) = C (ui,...,u;i+\Vui,...,Vu^+^) 

where 

c {zi ,... ,Zd+i,^i,... ,^d+i) = det 

is holomorphic (complex analytic). Any other constraint that can be written in 
terms of a holomorphic function f can be substituted to F. 


Zl ... Zd+l 1 \ 
Cl ■■■ Cd+ijy'’ 


In view of ([^ and (101, the following conclusion then naturally follows. 

Corollary 7. Assume that (§ and ( |11[ ) hold true. When d = 2 assume that 
n is convex and when d > 3 assume that holds. Suppose that A Q D and take 
D' (£ n. Then 


{(wi, ... ,u}d+i) € : {wfcjfc X {l,a;i... ,Xd} is F-complete in D'} 

is open and dense in . 


Remark 8. This result shows that almost any d + 1 frequencies in A give a F- 
complete set. An a priori estimate on the lower bound C for F (as in Theorem]^ 
cannot be obtained for arbitrary frequencies in an open and dense set in as 

this bound tends to zero when the frequencies are chosen near the residual set. 


It is easy to verify that this result is optimal in dimension 1 and 2 numerically. 
To intuitively see that AT = d + 1 is natural in any dimension, consider the level set 
uIj = 0, for a given lu. It is a priori a d — 1 dimensional object. If we consider the 
intersection of two such level sets for wi and uj 2 , we expect the resulting object to 
be d — 2 dimensional, the intersection of d such level sets to be zero dimensional, i.e. 
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Figure 1. A sketch of the d + 1 rationale in dimension 1 and 2. 
On the left, we represent two solutions with the same boundary 
data when d = 1: the locus {u^j =0} moves with the parameter 
uj, and the intersection of two zero sets is empty. On the right, we 
represent three zero level sets of the function F in two dimensions. 
The outer circle represents the boundary of O. When only two 
values of w are used, the zero level sets still contain common points, 
but the intersection of the three level sets is empty. 


discrete, and the intersection of (i+ 1 level sets to be empty. Figure [^is a graphical 
illustration of this idea. 

The rest of the paper is devoted to the proof of Theorem 


2. Proof of Theorem |4] 

First recall that the analyticity of the coefficients implies the analyticity of the 
solutions. Let C''^(r2; C) denote the space of complex-valued real analytic maps over 

n. 


Lemma 9 ([H]). Assume that (i and ( |11[ ) hold true. If lo G D and /* G 
H^dn-^R) then < e C'^(rf;C). 

Theorem 1^ will be a consequence of the following result. 

Proposition 10. Let 12,12' C be smooth domains such that 12'd 12. Let D C C 
be an open set such that 0 G D and AF D. Consider a map 

6»: C^(12;C), uj ^ 9^ 

such that for all x G LI, uj G D 9,^{x) G C is holomorphic and 9o(x) ^ 0. The set 

|(a;i,... ,Wd+i) G : imn(|6»<^J H-h \9u:a+i |) > o| 

is open and dense in . 

First, let us see why Theoremfollows from this result. 
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Proof of Theorem [7| Consider the map 

9: uj^ F =Mi,det 








In view of Lemma 9^^ G (7^(17; C). By the general fact that for any i and x 
the map lu G D (M^(a;), Vrt^(a;)) G is holomorphic [2, the map u! G D i-G 
9ui{x) G C is holomorphic for all x G 17. Moreover, 9q{x) 0 for all x G 17 by (12). 

We can apply Proposition and obtain that 

|(<Vi,... ,UJd+i) G : imn(|0^J + • • • + \9ui^^^ |) > o| 

is open and dense in 71'^+^. Note that the condition rninjy \9i^^ | > 0 is equiv¬ 
alent to 

for all X G 17' there exists k such that (wfc, /^,... \ > 0, 

which means that {oJk}k x {/*}i is a F-complete set of measurements in 17'. Indeed, 
defining 

f7'fc = {xGf7' : |F(ccfe,/\.../''+i)| >0}, 
we have 17' = Ufef7'^. This concludes the proof. □ 

The rest of this section is devoted to the proof of Proposition which is based 
on the structure of analytic varieties. 

An analytie variety in 17 is the set of common zeros of a finite collection of real 
analytic functions in 17, namely {x G 17 : 5 'i(x) = ••• = g^ix) = 0}, for some 
5i; ■ • ■ jffv G C''^(r7;C). For wi,... ,ujn G A we shall consider the analytic variety 

N 

Z{uJi,.. .,u}n) = {x G n : 6 »<^i(x) = • • • = 6»^„(x) = 0} = Q2'(a;j). 

Analytic varieties can be stratified into submanifolds of different dimensions. 

Lemma 11 f|22|h Let X be an analytic variety in 17. There exists a locally finite 
collection {A;}; of pairwise disjoint connected analytic submanifolds of LI (satisfying 
Whitney’s conditions) such that 

X = \JA,. 

i 

The decomposition X = U;A; is called a Whitney stratification of X. With this 
in mind, we can define the dimension of an analytic variety X = UjA; by 


(13) 


dimX := max dim A;. 


The main result leading to the proof of Proposition [T^ is the following 


Lemma 12. Under the hypotheses of Proposition 10 let LI” be a smooth domain 
such that LI” d 17 and X be an analytic variety in 17 such that X H 17" 7 ^ 0. Then 
the set 

{w G A : dim(Z(a;) n A n 17") = dim(A n 17")} 

is finite. 
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Proof. By contradiction, suppose that the set is infinite. Since A is compact, there 
exist G A, uJn ^ oj such that dim(Z(a;„) n X n ^2") = dim(X n ^2") and 

uJn ^ oj for all n G N. 


Therefore, in view of (131, for each n there exists a non-empty connected analytic 


submanifolds Sn such that 

(14) 

and 


Sn Q z{ojn) nxnn” 

dimS'ra = dim(X n fl”). 

for all n G N. Up to a subsequence, we have Xr, 


(15) 

Choose an arbitrary Xn G S, 
for some x G U". By Lemma 
U of X in n and a finite collection {Ai}i of analytic submanifolds of U such that 

xnu = UiAi. 
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X, 

applied to X, there exists an open neighborhood 


Moreover, since Xn G S„ and Xn —t x, up to a subsequence we have S'„ n t/ 0 for 
all n G N. As Sn Q X, up to a subsequence (and relabeling the collection Ai) we 
have for all n 


(16) 


SnCMI C Ai. 


Since by (15), dim(S'„ OU) = dim(X n U), and Ai C X C\U, this implies 
(17) dim(5'„ n C/) = dim Ai. 


In view of ( |l4| ) we have 9ui„{y) = 0 for all y G S'„. Therefore, by ( [I^ , ( p^ , |T7l 
Theorem 1.2] and 9^^ G we obtain SnfiU = Ai, whence 


Ai C Z{ijJn), n G N. 

As Xn G Ai for all n, we have a; G Ai. Thus, since Z{ujn) is closed, we infer that 
X G Z{uJn) for all n G N, namely 9^^{x) = 0 for all n G N. Since w 9i^{x) is 
holomorphic, this implies 9q{x) = 0, which contradicts the assumptions. □ 


We are now in a position to prove Proposition [T0| 

Proof of Proposition\10[ Since the map (wi,..., uJd+i) >—t minf^dfl,^^ |-|-- • |) 

is continuous, the set G = {(wi,... ,u}d+i) G : minj^d^^^J H-l-I^Ldd+J) > 0} 

is open. It remains to show that G is dense in . 

Take (wi,... ^Qd+i) G A^^^ and e > 0. Let U" be such that Pt' d U" (s U. We 
equip with the norm 

IKwi,... ,a;d+i)|| = maxlwfel. 

k 

We now want to construct an element (wi,... ,a;d-i-i) G G such that 
(18) II (^1’ ■ • ■ ’ ^d+i) — (wi, ■ • •, Wd-i-i) II < e. 

If dimZ(a;i) < d—1, set wi = wi; obviously we have |a;i—wi] < s. If dimZ(a;i) = d, 
by Lemma [T^ we obtain that the set 

{lo G A '■ dim{Z{uj) C U") = d} 
is finite. Therefore, we can choose wi G A such that 

dim(Z(a;i) n U") < d — 1 

and \uji — (hiI < £. Suppose now that we have constructed (^i,... , 0 ;^ such that 
\ujj — L 0 j\ < e for all j = Let us describe how to construct uJk+i- If 
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Z(uji,... ,uJk) n = 0, then it is enough to choose LOk+i = ^k+i- Otherwise, 
applying Lemma 12 with X = Z(uJi ,..., cuk), we obtain that the set 


{ui & A'. dim(Z'(a;) n ^(wi,... jWfc) n Vt") = dim(Z(a;i,... ,uik) H O")} 
is finite. Therefore, we can choose uJk+i € A such that 

dim(Z(tLii,..., w/c+i) n O") < dim(Z(wi,..., ujk) H Vt”) 

and \ijJk+i — Cjk+i\ < E. Therefore, as we have dim(Z(a;i) H fl") < d — 1, we obtain 
dim(Z(wi,..., Wd+i) n O") < 0, namely Z{uji,... ,0Jd+i) H O" = 0. In other words. 


(wi,..., oJd+i) G G. By construction, (181 is satisfied. This concludes the proof. □ 


3. Conclusions 


In this work we have showed that, under the assumption of real analytic coeffi¬ 
cients, almost any d -|- 1 frequencies in a fixed range give the required constraints, 
where d is the dimension of the ambient space. The proof is based on the structure 
of analytic varieties, and so the hypothesis of real analytic coefficients is crucial. 
To prove (or disprove) this result under weaker hypothesis on the coefficients a 
different approach is required. 

While this result seems optimal for an a priori fixed number of boundary con¬ 
ditions and for a somewhat arbitrary constraint function F, it could be that less 
than d -|- 1 frequencies are required if more boundary conditions are allowed (e.g., 
a set of d frequencies and d x (d -I- 1) boundary conditions to choose from). 
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